Charge imbalance in superconductors in the low-temperature limit 
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We explore charge imbalance in mesoscopic normal- metal/superconductor multiterminal struc- 
tures at very low temperatures. The investigated samples, fabricated by e-beam lithography and 
shadow evaporation, consist of a superconducting aluminum bar with several copper wires forming 
tunnel contacts at different distances from each other. We have measured in detail the local and 
non-local conductance of these structures as a function of the applied bias voltage V, the applied 
magnetic field B, the temperature T and the contact distance d. From these data the charge- 
imbalance relaxation length Aq* is derived. The bias-resolved measurements show a transition from 
dominant elastic scattering close to the energy gap to an inelastic two-stage relaxation at higher 
bias. We observe a strong suppression of charge imbalance with magnetic field, which can be directly 
linked to the pair-breaking parameter. In contrast, practically no temperature dependence of the 
charge-imbalance signal was observed below 0.5 K. These results are relevant for the investigation 
of other non-local effects such as crossed Andreev reflexion and spin diffusion. 



I. INTRODUCTION 

Non-equilibrium phenomena in superconductors have 
been investigated intensively since the 1970s. Both ex- 
perimental and theoretical investigations of charge imbal- 
ance have focused mostly on temperatures near the criti- 
cal temperature T c of the superconductor. In this regime, 
charge imbalance is easily accessible to experiments ba- 
sically due to the divergence of the signal towards T c , 
and excellent theoretical approximations are available^— 
Recently, the investigation of non-local transport proper- 
ties of superconductors has gained new impetus from two 
separate but loosely related fields. One is the investiga- 
tion of spin-dependent transport? 4 - and in particular spin 
diffusion and accumulation in the context of spintron- 
ics. The second is the investigation of coherent non-local 
effects such as crossed Andreev reflection^ which might 
be useful for quantum information processing. In the dif- 
fusive quasi-onedimensional structures typically used for 
such experiments, the magnitude of the signals due to 
these phenomena scale with the normal-state resistance 
of the superconductor over a characteristic length scale, 
which is given by the charge-imbalance relaxation length 
Aq* ~ 10 /imp the spin-diffusion length A s f ~ 1 /im^r- 
and the coherence length £ ~ 0.1 /mijiS respectively. 
Consequently, the signals due to charge imbalance are of- 
ten the largest, and make an unambiguous identification 
of the other phenomena difficult. In particular for the 
investigation of crossed Andreev reflection, experiments 
far below T c are necessary. Despite the vast amount of 
experimental and theoretical literature on charge imbal- 
ance, surprisingly little is known about the subject at 
very low temperatures ^ This is probably due to the 
fact that no simple theoretical models are available for 
this regime, and that the interpretation of the widely- 
used non-local resistance measurement scheme becomes 
increasingly difficult as temperature is lowered. 




FIG. 1. (Color online) SEM image of sample A illustrating 
the experimental scheme. Five copper (Cu) fingers are con- 
nected by tunnel contacts to an aluminum bar (Al). For one of 
the possible injector-detector pairs the bias and measurement 
scheme used for charge imbalance detection is shown. 



In this article, we report on a detailed investiga- 
tion of non-local conductance rather than resistance 
in superconductor/normal-metal hybrid structures at 
temperatures T <C T c . The samples consist of a 
quasi-onedimensional superconducting wire, with several 
normal-metal tunnel junctions attached to it. From these 
experiments, we deduce the charge-imbalance relaxation 
length Aq» as a function of bias voltage, temperature and 
magnetic field. The results allow a detailed comparison 
to theoretical predictions for different charge-imbalance 
relaxation mechanisms, and and an assessment of the 
possible impact on the interpretation of experiments on 
crossed Andreev reflection and spin diffusion. 



II. THEORY 

Charge imbalance (CI) can be described using two 
different theoretical frameworks, the quasiparticle, or 
two-fluid, approach , 1 ' 2 ' 12 ! 13 and quasiclassical Green's 
functions^ The relation of these two approaches has 
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been discussed extensively in the literature! 15 ' 16 We sim- 
ply note here that the Green's function method is more 
general, and can be applied, e.g., to inhomogeneous 
superconducting states and situations with strong pair 
breaking. We will nevertheless use the quasiparticle ap- 
proach, due to its conceptual (and computational) sim- 
plicity, and discuss its shortcomings where necessary. 

We consider a quasi-onedimensional superconductor of 
length L along the x axis, with several normal-metal elec- 
trodes attached via tunnel junctions. These electrodes 
will serve both to inject noncquilibrium quasiparticlcs 
into the superconductor, and to detect them. The quasi- 
particle energies E are given by 



E= Ve 2 + A 2 , 



(1) 



where e is the normal-state electron energy relative to 
the Fermi energy, and A is the pair potential. The nor- 
malized quasiparticle density of states is 



n(E) = e(E-A)^- 



(2) 



where is the Heaviside function. Charge imbalance is 
defined as 



Q* = 2N 



q(e)f(e)de, 



(3) 



where Nq is the density of states per spin at the Fermi 
energy in the normal state, q(e) — e/E is the effective 
quasiparticle charge in units of the elementary charge e, 
and /(e) is the quasiparticle distribution function. By 
using e rather than E as independent variable, we keep 
track of the electron-like (e > 0) and hole-like (e < 0) 
branch of the quasiparticle spectrum. In thermal equi- 
librium, /(e) is given by the Fermi distribution fc(E). 

It is apparent from ([3]) that Q* is non-zero only if the 
populations of the electron- and hole-like branches are 
unequal, i.e., /(e) — /(—e) ^ 0. Besides this charge-mode 
(transverse) non-equilibrium, there is also an energy- 
mode (longitudinal) non-equilibrium characterized by 
/(e) + /(-e) - 2/tCE) ^ 0. As will be shown below, 
the latter enters transport properties only indirectly via 
the self-consistency equation for the pair potential, 



1 + V 



l-2/(e) 



A 2 



de = 0, 



(4) 



where V is the pairing interaction. 

The electric current through a tunnel junction between 
a normal metal held at bias voltage V and a noncquilib- 
rium superconductor is given by the sum of the usual " Gi- 
aever" tunnel current It(V)^ and a bias-independent 
extra current Iq* due to CI^ 



I{V) = h{V) + I Q 



(5) 



Here, 



where Gn is the normal-state tunnel conductance, and 
the Fermi functions describe the electron occupation in 
the normal metal. The excess current is given by 



GnQ* 
2eNn 



(7) 



As mentioned above, the current Iq* depends on the 
nonequilibrium distribution /(e) only via the charge 
mode, whereas It(V) indirectly depends on the energy 
mode via the gap equation (|4]). 

In the vicinity of the critical temperature T c of the su- 
perconductor, the deviation of the distribution function 
/(e) from equilibrium is small and can be approximated 
by a Fermi function with a shift A/i of the chemical po- 
tential of the quasiparticles relative to the chemical po- 
tential of the Cooper pairsJ^ The most convenient mea- 
surement technique in this situation is the widely used 
non-local voltage detection scheme, in which the voltage 
between a normal-metal detector junction and the su- 
perconductor is measured. The voltage adjusts such that 
I(Vdet) = 0, i.e., the current Iq* is cancelled by the back- 
flow ^(Vdet)- In the regime where the chemical-potential 
model applies, this voltage is given by eVdct = A/i, i.e., 
it directly measures the single parameter that charac- 
terizes CI. Voltage detection, however, is less useful in 
the low-temperature regime that we are interested in for 
two reasons. First, at low temperature the chemical- 
potential model breaks down, and the detector voltage 
has no longer a simple physical meaning. Second, the 
non-linearity and temperature dependence of the tun- 
nel current It distort the measured signal. This has al- 
ready been noted in the earliest experiment on charge 
imbalance^ where the raw data at low temperature had 
to be corrected by the temperature dependence of the 
detector junction for comparison with theory. To avoid 
these shortcomings, we measure the detector current 
-fdet(Vdct = 0) = Iq*. Here, the only detector property 
that enters is the (constant) normal-state conductance, 
and therefore idet is a direct measure of Q* even for arbi- 
trary non-equilibrium distributions. Experimentally, we 
measure the differential non-local conductance, and we 
will now derive the expression used to evaluate our re- 
sults. 

The distribution function /(e) is driven out of equi- 
librium by tunnel injection into a fixed volume f2 of the 
superconductor at a rate given b y 12 ' 18 



dt 



1 f 1 



7tu 



-l^-^W-ME + eV))}, (8) 



where r tun = Gn /2N Qe 2 . Here, we use the full 
non-equilibrium distribution /(e) in the superconduc- 
tor, rather than the thermal distribution used in Ref. 
[l8l . However, we retain the assumption that the normal- 
metal electrode is at thermal equilibrium, described by 



It(V) = ^ J n(E) (f T (E - eV) - fc{E + eV)) dE 
o 

(6) 
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the shifted Fermi functions fc{E ± eV). It is custom- 
ary to define an injection efficiency F* — eilQ* /I. Since 
we will be interested in the differential conductance, we 
use the spectral quantity f*(E) = e 2 /E 2 rather than 
the usual integral definition. The expression for /* ne- 
glects the non-equilibrium contributions in ((5]) and ((8|), 
which are unimportant for an injector junction biased 
at \eV\ > A. On the other hand, a detector junction 
held at V — leads to charge imbalance relaxation at a 
rate r t ~^, and care must be taken to ensure that this 
rate is negligible compared to bulk relaxation mecha- 
nisms for non-invasive detection* 1 ^ Once injected, non- 
equilibrium quasiparticles diffuse along the wire with an 
energy-dependent diffusion constant D(E) = v g Dn , 16]20 
where £>n is the normal-state diffusion constant, and 
v s = \e\/E is the normalized group velocity of quasipar- 
ticles. Concomitantly, the non-equlibrium distribution 
relaxes due to different mechanisms, including inelastic 
electron-phonon scattering ) 2 ' 12 ! 21 elastic impurity scat- 
tering in the presence of gap anisotrop y 12 i 22 and mag- 
netic pair breaking j 14 i 23 ' 24 Charge- imbalance relaxes over 
a characteristic time tq* = Q*/Q*, which we will also 
assume to be energy-dependent, leading to an exponen- 
tial relaxation on the length scale Aq* = yj Dtq„ . The 
steady-state distribution is achieved when injection and 
relaxation rates are equal. Assuming that the effective 
injection volume is given by a wire section of length 2Aq* , 
we find that the non-local conductance due to charge im- 
balance is given by 

where Gi n j and Gdet are the normal-state conductances 
of the injector and detector junctions, is the normal- 
state resistivity of the superconductor, and A is the cross- 
section of the wire. The factor g* accounts for ther- 
mal smearing, injection efficiency, etc., and is of order 
unity. At T = 0, and neglecting energy relaxation, 
g* = n(E)f*/v s = 8(eVi nj - A). Equation © will form 
the basis for our data analysis. 

III. EXPERIMENT 

Fig. [T] shows the relevant part of one of the three sam- 
ples discussed, together with a scheme of the measure- 
ment setup. All samples are fabricated by e-beam lithog- 
raphy and shadow evaporation techniques. In the follow- 
ing, the processing details are explained with the help of 
the sample parameters listed in Table HI 

First, a copper film with thickness icui = 25 — 30 nm 
is evaporated onto a thermally oxidized silicon substrate. 
This first layer will form Ohmic interconnections to the 
subsequent layers. In a second evaporation step, the 
superconductor, an aluminum bar of thickness <ai and 
width wai, is deposited under a different angle, shift- 
ing the design to create intended overlaps only. In order 
to provide the formation of an insulating layer on the 



TABLE I. Characteristic parameters of the three samples A- 
C. Aluminum film thickness £ai, width wai, and normal state 
resistivity pAi at T = 4.2 K, range of contact distances d and 
contact conductances G, critical temperature T c , critical field 
B c and energy gap Aq. 





£ai 


WAI 


Pai 


d 


G 


T c 


Be 


A 




(nm) 


(nm) 




(ttm) 


OiS) 


(K) 


(T) 


(peV) 


A 


30 


140 


4.9 


1-12 


230-270 


1.39 


0.53 


208 


B 


30 


190 


4.9 


0.2-9.7 


260-350 


1.38 


0.58 


218 


C 


12.5 


140 


11.1 


0.5-6.5 


370-490 


1.5 


1.73 


225 



top, the aluminum is then oxidized in situ by applying 
the equivalent of 1 Pa of oxygen for 10 min. In the fi- 
nal evaporation step, a second layer of copper (£cu2 = 
30 nm) is deposited under a third angle forming five tun- 
nel contacts with the aluminum. The contact distances 
between neighbouring copper fingers presented in sample 
A are about 1 /im, 2 /xm, 4 /im and 5 /xm from left to 
right, respectively. 

For the transport experiment the samples are mounted 
into a shielded box thermally anchored to the mixing 
chamber of a dilution refrigerator. The measurement 
lines are fed through a series of filters to eliminate rf and 
microwave radiation from the shielded box. A voltage 
V cx consisting of a dc bias and a low-frequency ac excita- 
tion is applied to the injector contact, and the ac part of 
the resulting current L ln j is measured with a lock-in tech- 
nique. Simultaneously, the ac current T^et is measured 
through the second contact, the detector. The local and 
non-local differential conductances gi n j — dl- m j / dVi n j and 
g„i = d/dct / dVi n j are extracted from the ac signals. Volt- 
age and current polarities are indicated in Fig. Q] by plus 
signs and arrows, respectively. All contacts are measured 
in a three-point configuration with a series resistance of 
about 90 f2 coming from the measurement line. 

IV. RESULTS 

A. Contact and film characterization 

To characterize the tunnel contacts and the properties 
of the superconducting film, we first discuss the tunnel 
spectra of the individual junctions. Fig. [5] shows the lo- 
cal differential conductance <7i n j as a function of bias Vl n j 
for one of the tunnel junctions of sample B. Panel (a) 
represents the temperature dependence at zero magnetic 
field, whereas panel (b) depicts the variation with mag- 
netic field applied in the plane of the substrate along the 
Cu wires for constant temperature T = 50 mK. At lowest 
temperature and zero field, the differential conductance 
is completely suppressed at low bias, with sharp peaks 
at the energy gap, showing the high quality of the oxide 
tunnel barrier. Upon increasing the temperature or the 
magnetic field, the features are broadened and the gap is 
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FIG. 2. (Color online) Local differential conductance g — 
dl /dV of one contact of sample B as a function of bias voltage 
V for (a) different temperatures T and (b) different applied 
magnetic fields B. Symbols represent measured data, lines 
are fits to the model described in the text. 



reduced. Since we are particularly interested in the de- 
pendence on magnetic field, we have used a slightly more 
elaborate model than (|6|) to fit the data. In the presence 
of a magnetic field, the spin-resolved density of states in 
the superconductor can be described by2£ 



n±{E) = iRe 




(10) 



where the complex quantities u± have to be determined 
from the implicit equation 



u± 1 



r 



i 



i 



Here, /zb is the Bohr magneton, T is the pair-breaking pa- 
rameter, b so — h/3r so A measures the spin-orbit scatter- 
ing strength, and we have dropped a small higher-order 
term. The fits in Fig. [5] were obtained by replacing the 
BCS density of states n{E) by n + (E) + n^(E) in @. 
During fitting, T and B were taken from the experiment, 
and the remaining parameters were varied. Including the 
Zeeman splitting was found to be necessary for the data 
at higher fields. The small spin-orbit term gave minor 
improvements of the fits, but could not be determined 
precisely. We simply chose a suitable value b so ~ 0(0.01) 
at high field, and kept it fixed for all other fits. Similar- 
values can be found in the literature^ The quality of the 
fits was excellent for samples A and B, as shown for the 
latter in Fig. [21 For the thin- film sample C, the quality 
of the fits was rather poor, and no reliable parameters 
could be obtained. 

For samples A and B, where the fits were reliable, 
we show the pair-breaking parameter T as well as the 
pair potential A as a function of the magnetic field in 
Fig. [U In panel (a), the normalized pair potential A/Ao 
is displayed as a function of B/B c , together with the ex- 
pectation ln(A/A ) = -(7r/4)(r/A)j21 The critical pair- 
breaking strength for the suppression of superconductiv- 
ity is given by 2T = Apr - and together with T oc B 2 for 



FIG. 3. (Color online) (a) Pair potential A as a function of 
the applied magnetic field B. A is normalized to its zero- 
field value Ao, and B is normalized to the critical field B c . 
(b) Normalized pair- breaking parameter T/Ao as a function 
of the (B/B c ) 2 . The lines are predictions from pair-breaking 
theory. 
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FIG. 4. (Color online) a) Differential conductance gdet of the 
left-most contact of sample A as a function of bias voltage 
Vdot, with additional voltage bias Vi n j applied to the neighbour 
contact at 1 /im distance, b) Normalized pair potential A/Ao 
as a function of injector bias Vinj . The line is a guide to the 
eye. 



a thin film in parallel magnetic field,— we can rewrite T 
as 



r 



1JB 

2 \B r 



(11) 



Fig. [3(b) shows r/A as a function of (B/B c ) 2 . The 
solid line is the theoretical expectation (fTTj) . As can be 
seen, samples A and B can be described perfectly well by 
standard pair-breaking theory, and we will assume (jllj) 
to hold also for sample C later on. 



B. Energy-mode non-equilibrium 

While we are mostly interested in charge imbalance, 
we have also investigated the impact of energy- mode non- 
equilibrium in our samples. To this effect, we monitor the 
differential conductance <Jdet of a detector junction while 
non-equilibrium quasiparticles are injected through a sec- 
ond nearby injector contact. Fig. 0(a) shows the differen- 
tial conductance g^et of the left-most contact of sample 
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FIG. 5. (Color online) Non-local differential conductance 
g n i — d/dot/rfVinj for an injector/detector pair of sample A as 
a function of bias voltage Vinj (a) for different temperatures 
T and (b) for different applied magnetic fields B. 

A (see Fig. []} as a function of the local bias voltage Vdet 
for different injector bias Vi n j applied to the neighbouring 
contact at a distance of 1 fxm. We focus here only on the 
bias region of the gap features. As the injector bias is 
increased, the density-of-states peak in the conductance 
shifts to lower bias, and broadens slightly. The increased 
broadening of the gap features might signify an increased 
temperature of the normal-metal side of the junction due 
to the injection of non-equlibrium quasiparticles tunnel- 
ing out of the superconductor^ An alternative expla- 
nation would be an increased lifetime broadening of the 
density of states of the superconductor due to scattering 
of non-equilibrium quasiparticles. From our data, we can 
not make a clear decision between these scenarios. The 
evolution of the pair potential A as a function of injec- 
tor bias is plotted in panel (b), normalized to its value 
Ao at Vinj = 0. No significant change is observed for 
eVinj < A . As soon as eVi n j exceeds A « 200 /xeV, 
the energy gap drops quickly by a few percent, and con- 
tinues to decrease more slowly for higher bias. The gap 
reduction can be understood from the inspection of the 
self-consistency equation (fj|. For injector voltages in the 
vicinity of A, a large number of quasiparticles are in- 
jected due to the divergence of the density of states. In 
addition, these quasiparticles are very efficient in reduc- 
ing the gap due to the energy denominator in the inte- 
gral. Therefore, the initial decrease is steep, and then 
becomes more shallow as the density of states flattens, 
and the denominator increases. 



C. Charge-mode non-equilibrium 

Figure \5\ displays the non-local conductance g n \ as a 
function of the injector bias T4nj for one injector/detector- 
pair of sample A. Panel (a) shows data for different 
temperatures T without an applied magnetic field. At 
T = 50 mK, the non-local conductance is zero within the 
experimental resolution for bias voltages below A/e w 
200 /iV. Above the energy gap, the signal increases con- 
tinuously from zero with a finite initial slope up to a 



FIG. 6. (Color online) (a) Normalized nonlocal differential 
conductance gfni/GinjGdet as a function of injector bias voltage 
Vinj for different contact distances d. (b) Semi-logarithmic 
plot of gni/GinjGdot as a function of contact distance d for 
different injector bias Vi n y The solid lines are fits to ((9|. 

broad maximum at Vinj ~ 450 /iV before it decreases 
slowly again. With increasing temperature, the signal 
smears out around the gap, whereas the value at high 
bias remains unchanged. Panel (b) shows the impact of 
a magnetic field B applied in the substrate plane along 
the direction of the copper wires. In contrast to temper- 
ature, the signal depends strongly on the magnetic field. 
The initial slope decreases, and the maximum decreases 
and shifts to higher bias until it is no longer observable 
within our bias range for B > 100 mT. 

Figure [SJa) shows the non-local differential conduc- 
tance for several injector/detector contact pairs of sam- 
ple A at T = 50 mK and B = 0. Since g n \ oc Gi n jGd e t, 
we have normalized the data accordingly to exclude the 
impact of small variations of the junction conductances. 
The overall signal magnitude decreases with increasing 
contact distance while the shape remains unchanged. 
Panel (b) shows the normalized non-local conductance 
as a function of contact distance d for different injector 
bias on a semi-logarithmic scale. The solid lines are fits 
to the exponential decay predicted by . The quality of 
the fits is generally good, except for the very small sig- 
nals at lowest bias. From the fits, the relaxation length 
Aq* and the amplitude 



can be extracted. 

The charge imbalance relaxation length Aq* extracted 
from these fits is shown in Fig. [7] as a function of injector 
bias for different temperatures T (a) and magnetic fields 
B (b). The data resemble those of the non-local con- 
ductance shown in Fig. [5l This is not surprising, since 
the signal amplitude a is itself proportional to Aq* . A 
noticeable difference is that Aq* is nearly independent of 
temperature. This indicates that the temperature depen- 
dence seen in Fig. Oa) is mostly due to thermal broaden- 
ing of the distribution in the injector contact rather than 
a change in relaxation rates. In contrast, the suppression 
of the non-local conductance upon increasing the mag- 
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FIG. 7. (Color online) Charge imbalance relaxation length 
Aq* as a function of injector bias voltage Vi n j for (a) different 
temperatures T and (b) different applied magnetic fields B. 
The line is the result of a numerical simulation described in 
section W\ 



TABLE II. Maximum values of the relaxation length Aq* , 
the relaxation time tq* , and the ratio tq* /rtun for all three 
samples. 



sample 


Aq. 




TQ'/r tun 




(/xm) 


(ns) 




A 


5.2 


7.8 


0.01 


B 


4.3 


5.2 


0.01 


C 


3.0 


5.9 


0.05 



netic field is reflected in the pronounced field dependence 
of Aq. , indicating an increase of the relaxation rate. 

Similar results to those presented in Figs. [5]|7] were ob- 
tained for all three samples. From Xq* , we can calculate 
tq* = Aq, / Dj<jVg. The maximum values of Xq* and tq* 
obtained at lowest temperature and zero magnetic field 
are listed in table [TTl along with the maximum of the 
ratio tq* /Tt un . We find tq* / Tt U n -C 1 for all samples, 
confirming that our detector junctions are non-invasive. 

We will now focus in more detail on the suppression 
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FIG. 8. (Color online) (a) Normalized charge-imbalance re- 
laxation rate h/ Aotq* as a function of the normalized pair- 
breaking parameter T/Ao. The line is a guide to the eye. 
(b) Normalized charge-imbalance relaxation rate Ti/Ttq* as a 
function of normalized injector bias eVi n j/Ao. The line is a 
joint fit to the data of all three samples. 



FIG. 9. (Color online) g* as a function of normalized injector 
bias eVi n j/Ao (a) for all samples at B = and (b) for sample 
A at different magnetic fields B. 



of charge imbalance as a function of magnetic field. Fig- 
ure [5Ja) shows the normalized charge-imbalance relax- 
ation rate U/ Aqtq* as a function of the normalized pair- 
breaking parameter r/A for fixed injector bias. Here we 
have made use of equation to calculate V from B for 
all three samples. The data from all samples fall onto a 
single line, and the relaxation rate at zero field is negligi- 
ble on the scale of the plot. We note that in the magnetic- 
field range of the plot (B < 0.5B C ) the spectral properties 
of the superconductor remain almost unchanged. The re- 
duction of A, for example, is less than 10% in this range. 
The relaxation rate due to elastic pair-breaking perturba- 
tions such as magnetic impuritieS ) 14 ' 23 i 30 supercurrem 
and applied magnetic field 3 -^ has been calculated both 
within the quasiparticle description used by us, and from 
quasi-classical Green's functions. We note that by con- 
vention, rates from the Green's function formalism differ 
from those of the quasiparticle description by the factor 
When properly adjusted, the rate is predicted to 

be 
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1 



r A 2 
< . 

h Ee 



(13) 



where a is a numerical prefactor of order unity which we 
will use as a fit parameter. From linear fits of the data 
in panel (a) we can extract H/Ttq* . The result extracted 
from such fits is plotted in panel (b) as a function of 
normalized injector bias eV; n j/Ao for all three samples. 
The solid line is a joint fit of all data to (fTBf (where we 
have set E = eV; n j and A = A ), with a = 0.73. 

From the signal amplitude a extracted from the fits in 
Fig. [6j we can calculate the prefactor g* using (| 1 2[) to- 
gether with known sample parameters and Aq* extracted 
from the same fits. The results are plotted in Fig. [H] as a 
function of normalized injector bias (a) for all samples at 
B = and (b) for sample A at different magnetic fields. 
At B = 0, g* follows the expectation g* s» 9(eVi n j — A) 
for samples A and B, whereas it deviates at low bias both 
for sample C, and in the presence of a magnetic field. 
Since g* depends on a combination of several spectral 
properties of the superconductor, we can not identify the 
precise cause of these deviations. For sample C, the en- 
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hanced energy-mode non-equilibrium due to the reduced 
film thickness may play a role. However, in all cases 
g* w 1 at high bias, which justifies our choice of using a 
wire section of length 2Xq* as injection volume. 

V. DISCUSSION 

The bias dependence of the non-local conductance can 
be understood as follows: Charge-imbalance relaxation 
takes place mostly at energies close to A, since the coher- 
ence factors for scattering between the electron- and hole- 
like branches, both elastic and inelastic, vanish at higher 
energies. In addition, in the low-temperature regime the 
inelastic contribution is expected to be negligible. 32 The 
elastic relaxation rate, both due to gap anisotropy and 
magnetic pair breaking, diverges for E — > A, and quickly 
drops at higher energies. Consequently, charge imbal- 
ance rises continuously from zero as the bias is increased 
above A, and in this regime relaxation is mainly due to 
direct scattering between the branches. As the bias is in- 
creased further, the direct relaxation rate decreases. On 
the other hand, energy relaxation due to inelastic scat- 
tering becomes important, and charge relaxation turns 
into a two-stage process^ Quasiparticles are first cooled 
to the vicinity of A by inelastic scattering, and are then 
scattered elastically between branches. If we assume that 
inelastic scattering is described by electron-phonon scat- 
tering in the Debye approximation, the inelastic rate 
quickly increases as more phonons become available at 
higher energy. Consequently, the non-local differential 
conductance begins to drop again after its initial increase. 
The bias at which the maximum appears scales roughly 
with the transition between the direct and two-stage re- 
laxation regimes. The temperature dependence, or lack 
thereof, can be understood easily from this picture. Elas- 
tic relaxation depends only weakly on temperature via A. 
For inelastic relaxation, at low temperatures the Bose 
factors for emission and absorption of phonons become 
« 1 and « 0, respectively. Relaxation is then dominated 
by phonon emission at a rate which only depends on 
the bias-dependent energy of the quasiparticles, but not 
on temperature. Consequently, the relaxation rate, and 
Xq* , become practically independent of temperature over 
the entire bias range, as observed in Fig. Ufa). We have 
identified the quasiparticle energy with the bias voltage 
throughout our data analysis in section IPV1 This has to 
be taken with some caution due to the presence of inelas- 
tic energy relaxation. Our approximation mainly affects 
the calculation of tq* = Xq, / 'D^v g , where v g should ac- 
tually be an average over energy. However, since inelastic 
relaxation is important mostly at higher energies, where 
v g ss 1, we assume that the error is small. This is corrob- 
orated by the observation that the position of the max- 
imum in g n \ does not depend much on contact distance, 
as seen in Fig. [3a). We would nevertheless like to stress 
that Xq* and tq* are not really spectral quantities, but 
depend in detail on the non-equilibrium distribution and 



bias conditions. 

To make a quantitative connection to microscopic the- 
ory, we have performed numerical simulations of the 
quasiparticle Boltzmann equation, basically following 
Ref. |22|. We have used the onedimensional form of the 
Boltzmann equatio n 16 ' 22 

df d 2 f 

~fo~ v s DN d^ = 7tun - 7ci - Tin- (14) 

The elastic and inelastic relaxation rates 7 e i and 7; n are 
given by equations (2.16) and (2.17) of Ref. respec- 
tively. The injection rate is given by ©. Steady-state 
solutions /(e, x) of the Boltzmann equation (|14l) were ob- 
tained by numerical iteration on a discretized grid. Both 
injector and detector junction were included on an equal 
footing, with injection rates and currents given by (|8]l 
and (|5|). Here, the injection volume £1 is given by the 
grid point size, and the spreading of charge imbalance 
over the length scale Xq* is included microscopically in 
the diffusion term. The granularity of the grid was cho- 
sen sufficiently small (Ae = 20 /J.eV, Ax = 500 nm) not 
to affect the results. Parameters such as wire geometry, 
diffusion constant, contact conductances Gj n j and Gdetj 
etc., were taken directly from the experiment, leaving 
only the characteristic electron-phonon scattering time 
To and the average gap anisotropy (» 2 ) as free parame- 
ters. The results of the simulation are shown as a solid 
line in Fig. EJb), where we have inserted the typical val- 
ues r = 100 ns and (a 2 ) Q = 0.03 from Ref. 22. The 
overall magnitude and shape is predicted correctly by 
the simulation, with Xq* ps 5 /im, and a broad maxi- 
mum at Vi n j fa 400 /xV. However, in detail the agreement 
is rather poor. The slope is too steep, both below and 
above the maximum. The relaxation rate due to impu- 
rity scattering has been calculated both within the quasi- 
particle approach ) 12 i 22 and using quasiclassical Green's 
functions^ 2 - In contrast to the quasiparticle approach, the 
Green's function approach predicts different energy de- 
pendences for the clean and dirty limits. This difference 
might explain the discrepancy at low bias. At high bias, 
electron-phonon scattering in the Debye approximation 21 
apparently overestimates the energy dependence of the 
relaxation rate. A weaker energy dependence has been 
predicted for electron-phonon scattering in the presence 
of disorder i^ 3 . - — Also, it has been argued that disorder- 
enhanced electron-electron interaction may dominate in 
particular in aluminum with its weak electron-phonon 
interaction^ 6 - - — A simulation based on microscopic the- 
ory for both electron-electron and electron-phonon inter- 
action might provide detailed insight into the relaxation 
mechanisms here, but this is beyond the scope of this 
article. 

We now focus on the dependence on magnetic field. 
Magnetic pair breaking adds an elastic contribution to 
the relaxation rate i 14 ' 23 ' 24 ' 30 ' 31 Consequently, as B in- 
creases, the initial slope of the differential conductance 
decreases, and the maximum, i.e., the transition to the 
two-stage relaxation regime, shifts to higher bias. We 
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find that the magnetic contribution to the relaxation rate 
is directly proportional to the pair-breaking parameter 
r, and that its energy dependence follows the theoret- 
ical prediction independent of sample details. Our ob- 
servation Tq} oc r is markedly different from the well- 
established approximation Tq} oc \/l + 2teT/H valid for 
T — > T c ^ where Tg is the energy relaxation time (note 
that our T/h is the mag netic pair-breaking rate, i.e., the 
quantity t~ 1 of Ref . Il4l ) . Magnetic relaxation dominates 
all other contributions even at very low magnetic fields, 
where the spectral properties of the superconductor are 
almost unaffected by pair breaking. 

We finally discuss the possible impact of charge imbal- 
ance on the observation of other phenomena. We first 
note that the non-local conductance at subgap energies 
remains negligible at lowest temperatures, even with an 
applied magnetic field. This is not surprising, since even 
in the presence of magnetic pair breaking the supercon- 
ductor still has a well-defined energy gap for quasipar- 
ticle excitations, at least as long as pair breaking is not 
too strong. Also, the relaxation length is always larger 
than the coherence length (£ 100 nm for our sam- 
ples). Thus, the dependence of non-local conductance 
on bias or contact distance remains a good criterion to 
distinguish coherent subgap transport from charge im- 
balance. On the other hand, the observation of spin- 
dependent quasiparticle transport necessarily involves in- 
jection at energies above the gap, and ferromagnetic elec- 
trodes must be used. At magnetic fields of ~ 100 mT, 
which are easily reached by the fringing fields of elec- 
trodes made of elementary ferromagnets, Xq* can already 



be as small as 1 /im. This is similar to the spin-diffusion 
length A s f in aluminium^"* and consequently great care 
must be taken to distinguish charge imbalance and spin- 
dependent transport. 



VI. CONCLUSION 

In conclusion, we have presented a detailed investiga- 
tion of charge imbalance in superconductors in the low- 
temperature regime. From our measurements, we have 
extracted the charge-imbalance relaxation length as a 
function of bias voltage, temperature, and magnetic field. 
The bias- dependent results allow for a detailed compari- 
son with different relaxation mechanisms. In particular, 
we have shown a transition from dominant elastic relax- 
ation in the vicinity of the energy gap to an inelastic 
two-stage relaxation at high bias. The dependence on 
magnetic field follows theory with remarkable accuracy, 
and is clearly different from the known approximations 
for the high-temperature regime. The strong reduction 
of the relaxation length with magnetic field has possi- 
ble implications for the interpretation of spin-diffusion 
experiments. 

During the preparation of this manuscript we became 
aware of three related studies of charge imbalance at 
low temperatures^r— We thank K. Yu. Arutyunov, 
W. Belzig and D. S. Golubev for useful discussions. This 
work was supported by the DFG Center for Functional 
Nanostructures, and by the DAAD within the RISE pro- 
gram. 
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